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Abstract. The presence of a white dwarf in a resolved binary system, such as Sirius, 
provides an opportunity to combine dynamical information about the masses, from as- 
trometry and spectroscopy, with a gravitational red-shift measurement and spectropho¬ 
tometry of the white dwarf atmosphere to provide a test of theoretical mass-radius rela¬ 
tions of unprecedented accuracy. We demonstrated this with the first Balmer line spec¬ 
trum of Sirius B to be obtained free of contamination from the primary, with STIS on 
HST. However, we also found an unexplained discrepancy between the spectroscopic 
and gravitational red-shift mass determinations. With the recovery of STIS, we have 
been able to revisit our observations of Sirius B with an improved observation strat¬ 
egy designed to reduce systematic errors on the gravitational red-shift measurement. 
We provide a preliminary report on the refined precision of the Sirius B mass-radius 
measurements and the extension of this technique to a larger sample of white dwarfs in 
resolved binaries. Together these data can provide accurate mass and radius determina¬ 
tions capable of testing the theoretical mass-radius relation and distinguishing between 
possible structural models. 


1. Introduction 

Evolutionary model calculations (e.g. Wood et al. 1995; Fontaine et al. 2001; Berg¬ 
eron et al. 2001; Althaus & Benvenuto 1998, 1997) have established that the theoretical 
radii for DA white dwarfs are systematically dependent on stellar temperature, H layer 
mass and core composition. However, basic tests based on analyses of a combination 
of Baimer-line measurements and Hipparcos parallaxes (Vauclair et al. 1997; Barstow 
et al. 1997; Provencal et al. 1998) are restricted to the narrow range of masses of the 
isolated white dwarfs (near 0.6 M 0 ) and to the cooler end of the DA white dwarf T e n 
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range (<20,000K), where the M-R relation is less sensitive to temperature/layer mass 
and core effects. The hotter stars in the Hipparcos sample (e.g. Feige 24, G191-B2B), 
which lie at the greatest distances, have parallax uncertainties too large to provide a 
really thorough examination of the effects of differences in the core and envelope com¬ 
positions. ffowever, the presence of a bright companion, where the white dwarf resides 
in a binary system, yields a more accurate parallax, ffence, Sirius-like binary systems 
present a good opportunity of studying the M-R relation for a range of white dwarf 
temperatures and masses. 

In HST observing cycle 12 we obtained a high-S/N optical spectrum of the Balmer 
lines of Sirius B with HST/STIS, to provide precise measurements (the first with mod¬ 
ern electronic detectors) of T e ff, log g and gravitational redshift and test the theoretical 
white dwarf mass-radius (M-R) relation (Barstow et al. 2005). In cycle 19 we were 
granted more time to extend this work. We observed Sirius B again, to better under¬ 
stand potential systematic errors in the measurement of its gravitational redshift and 
also studied five other white dwarfs in resolved Sirius-like systems. We have obtained 
the first precise measurements of T e ff, log g and radial velocity for these WDs and 
combined them with the distances (from Hipparcos) and the primary radial velocities 
to measure the WD gravitational redshift. As a result we have obtained the stellar mass 
and radius, to test the M-R relation at high temperatures across a range of white dwarf 
masses. 


2. Observations 

We have followed the basic instrumental setup used in Barstow et al. (2005), and taken 
observations using the G430L and G750M gratings to cover the full Balmer line series. 
With its higher resolution, the G750M grating also provides an accurate radial velocity 
measurement (using the Her lines) and, as a result, gravitational redshift. In the original 
programme a 52x0.2" slit was used but this appeared to generate systematic uncer¬ 
tainties in both knowledge of the flux (possible light losses from the slit) and velocity 
measurement (uncertainty in positioning of the stellar image in the slit). Therefore, we 
used a wider slit for the G430L observation to ensure an accurate flux measurement and 
recorded a pair of observations in the G750M grating using a wide slit, 52x2.0" for the 
photometry and normalisation of the data, and a narrow slit, 52x0.05" to obtain a more 
accurate radial velocity. 

We observed a total of five white dwarfs in Sirius-like systems, SiriusB, HD2133B, 
14 Aur BCB, HR1358B and HD217411B. HD217411 was found to be a hierarchical 
triple system, with the white dwarf lying in an unresolved binary companion to the 
main sequence star. Therefore, an uncontaminated Balmer line spectrum could not 
be obtained for the white dwarf and no M-R measurement made. Our findings are 
published in (Holberg et al. 2014). 

All of the spectra were reduced with the calstis pipeline in iraf and analysed as 
described in (Barstow et al. 2005). Briefly, we used the spectral fitting code xspec in 
conjunction with a grid of pure hydrogen synthetic spectra, generated using tlusty and 
synspec, to obtain measurements of Teff and log g for each white dwarf. The stellar 
radial velocity was obtained from the redshift of the model spectrum compared to the 
rest wavelength of H a. 
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Figure 1. The points plotted on the left are the masses derived from the log g, 
while the points on the right are the masses derived from the gravitational redshift. 
Both sets of radii are derived from the normalisation of the spectra. The squares are 
the normalisation from the G750M grating and the circles from the G430L normal¬ 
isation. The solid black line is the CO core thick H layer and the dotted line is the 
thin H layer at 25 000 K. 


3. Results 


There are three different approaches that can be taken to measure the mass and radius 
for a white dwarf, depending on the data available. The first involves simply taking 
the T e ff and log g as measured from the Balmer line spectrum obtained in the G430L 
grating and using them as an input to a theoretical mass-radius relation, such as those 
provided by Fontaine et at. (2001) to calculate the mass and radius of the white dwarf. 

The second measurement involves calculating the mass directly from the log g 
measurement as, 


GM 

8 ~ ~W 


( 1 ) 


where G is the universal gravitational constant, and the radius can be calculated from 
the normalisation to the model (as STIS spectra are accurately flux calibrated), and the 
distance to the system. This measurement is dominated by errors on the system parallax 
which is known to high accuracy for Sirius B, but much less well for the other systems 
in our programme. The ESA Gaia mission should produce new measurements before 
2017 which will improve the accuracy of our results using this method. 

The third method involves determining the gravitational redshift. This is done by 
measuring the radial velocity of the white dwarf from the H a line and subtracting the K 
and y velocities of the system. The gravitational redshift can then be used to determine 
the mass of the system as, 


V g r = 


0.636 M 


R 


( 2 ) 


The radius is again calculated from the normalisation of the model to the data. As we 
have two sets of data - the Balmer lines and the ff a line, we have two measurements of 
the radius. 

Figure 1 shows the mass and radius derived from the gravitational redshift (right 
hand error bars) and from the log g (left error bars), both combined with the photomet- 
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Mass (M 0 ) 

Figure 2. The best-fit mass and radius (black dot with 1 and 2 cr error ellipses) 
compared with the astrometric constraints (box) and theoretical M-R relationships 
for a thick H layer (solid curve) and thin H later (dashed curve). The input parameters 
values are in the top right corner of the figure. 


ric normalisation. The square symbols are where the normalisation is from the G750M 
grating and the circles are from the G430L grating. The diagonal black lines are theo¬ 
retical mass-radius relations for CO core models at 25 000 K with thick (solid) and thin 
(dashed) H layers. It can be seen, as also found by (Barstow et al. 2005), that these two 
methods of deriving the masses do not agree with each other. Therefore, the refined 
observational approach designed to remove systematic errors has not resolved the issue 
and that some problem remains. 

Nevertheless, while there is still an analysis issue to be resolved, it is possible to 
combine all the methods of determining mass and radius and use a chi-squared minimi¬ 
sation to find the best estimate of these parameters. This is illustrated in figure 2, which 
shows the best fit (black dot) and the 1 and 2 cr error ellipses. The diagonal curves 
are as described in figure 1 and the box defines the astrometric limits (Bond, private 
communication). The input parameter values are in the top right corner of the figure. 
This illustrates that the precision of our measurements allows us to distinguish between 
theoretical M-R relations assuming different H layer masses. 

Figure 3 shows the mass-radius relation for 14 Aur C, Sirius B and HD2133 mea¬ 
sured as part of this programme. The measurements for Sirius B are derived from T e ff 
and log g and combined with a mass-radius relation from Fontaine et al. (2001). The 
mass has increased from our previous work (Barstow et al. 2005), but it can be seen that 
our eiTors have decreased significantly. The measurements for the other stars are from 
the log g derived mass, and the errors on the parallax dominate. At this stage, we do 
not have all the information for the other stars that we have for Sirius B. However, we 
can illustrate the potential outcome of the programme using a theoretical mass-radius 
relation from (Fontaine et al. 2001) to anchor the mass and radius measurements. Fig¬ 
ure 3 shows the results for 14 Aur C, Sirius B and HD2133, compared to curves for a 
range of temperatures and thick and thin H layer masses. For HD2133 and 14Aur-BC, 
the eiTor bars are dominated by the current uncertainties in the system parallax. Never¬ 
theless, the HD2133 result is capable of distinguishing between the models, helped by 
the increased dynamic range of the models at lower masses. Two results are plotted for 
Sirius B, from (Barstow et al. 2005) (left) and this work (right). The mass estimate has 
increased from our previous work, but it can be seen that our errors have also decreased 
significantly, allowing clear constraints to be placed on the models. 
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Figure 3. Masses and Radii derived from our STIS data. The Fontaine models 
(Fontaine et al. 2001) for CO cores and thick and thin H layers are shown for a range 
of T e ff. The result for Sirius B in Barstow et al. (2005) is also shown to the left of 
our new point which has smaller errors. 


4. Conclusions 

We have presented preliminary results from our F1ST/STIS observations of a sample 
of hot DA white dwarfs residing in resolved Sirius-type binary systems, with the aim 
of testing the theoretical mass-radius relation under a range of conditions. We have 
refined our observational technique, but this has not eliminated disagreements in the 
measurements of M and R from different approaches. It seems most likely that the 
issue is in the stellar atmosphere models used to fit the B aimer lines rather than the 
observations themselves. We have demonstrate the ultimate accuracy of our approach 
but further work is needed to refine some observations, such as the system parallaxes, 
and obtain more dynamical data for these systems. 
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